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Testing nonlinear vacuum electrodynamics with Michelson interferometry 
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We discuss the theoretical foundations for testing nonlinear vacuum electrodynamics with Michel¬ 
son interferometry. Apart from some nondegeneracy conditions to be imposed, our discussion applies 
to all nonlinear electrodynamical theories of the Plebariski class, i.e., to all Lagrangians that depend 
only on the two Lorentz-invariant scalars quadratic in the field strength. The main idea of the 
experiment proposed here is to use the fact that, according to nonlinear electrodynamics, the phase 
velocity of light should depend on the strength and on the direction of an electromagnetic back¬ 
ground field. There are two possible experimental setups for testing this prediction with Michelson 
interferometry. The first possibility is to apply a strong electromagnetic field to the beam in one 
arm of the interferometer and to compare the situation where the field is switched on with the 
situation where it is switched off. The second possibility is to place the whole interferometer in a 
strong electromagnetic field and to rotate it. If an electromagnetic field is placed in one arm, the 
interferometer could have the size of a gravitational wave detector, i.e., an arm length of several 
hundred meters. If the whole interferometer is placed in an electromagnetic field, one would have 
to do the experiment with a tabletop interferometer. As an alternative to a traditional Michelson 
interferometer, one could use a pair of optical resonators that are not bigger than a few centimeters. 
Then the whole apparatus would be placed in the background field and one would either compare 
the situation where the field is switched on with the situation where it is switched off or one would 
rotate the apparatus with the field kept switched on. We derive the theoretical foundations for these 
types of experiments, in the context of an unspecified nonlinear electrodynamics of the Plebariski 
class, and we discuss their feasibility. A null result of the experiment would place bounds on the 
parameters of the theory. We specify the general results to some particular theories of the Plebariski 
class; in particular, we give numerical estimates for Born, Born-Infeld and Heisenberg-Euler theories. 

PACS numbers: 03.50.Kk,11.10.Lm 


I. INTRODUCTION 

Modifying an earlier idea by Born PQ, in 1934 Born 
and Infeld ,2] suggested a nonlinear modification of vac¬ 
uum electrodynamics in order to get rid of the infi¬ 
nite self-energies of point particles that occur in the 
standard Maxwell theory. Their theory can be derived 
from a Lorentz-invariant Lagrangian. A few years later, 
Heisenberg and Euler [3j derived an effective Lagrangian, 
again Lorentz-invariant, from quantum electrodynamics. 
These are the two best known examples within the class 
of all Lorentz-invariant nonlinear electrodynamical theo¬ 
ries. More generally, Plebariski [4] and also Boillat [5] 
studied the whole class of nonlinear electrodynamical 
theories that can be derived from a Lagrangian depend¬ 
ing only on the two Lorentz-invariant scalars that are 
quadratic in the held strength. This class is often referred 
to as Plebariski nonlinear electrodynamics. For a review 
of the Born-Infeld theory we refer, e.g., to Bialynicki- 
Birula j6j. 
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The physical relevance of these nonlinear vacuum elec¬ 
trodynamical theories is being widely discussed in the 
literature. It is believed that at a certain held strength 
the Heisenberg-Euler deviations from standard Maxwell 
theory should be observable, and the Born-Infeld the¬ 
ory has gained increasing attention since it was realized 
by Tseytlin [7 that the Born-Infeld Lagrangian can be 
derived as an effective Lagrangian from some versions 
of string theory. Observable effects of (nonlinear) mod¬ 
ifications of the vacuum Maxwell equations have been 
discussed for many years, at least since the Ph.D. the¬ 
sis of Toll [5]. Up to now, the only effect predicted by 
such modified theories that has already been observed is 
light-by-light scattering (see 0); further experiments are 
under way, e.g. with the Large Hadron Collider at CERN 
m- There is also an ongoing experiment m aiming at 
verifying the birefringence in vacuo as predicted by the 
Heisenberg-Euler theory. Also, it might be possible to 
measure the influence of background fields on the prop¬ 
agation speed of light in the laboratory. For the case of 
the Born-Infeld theory, such experiments have been sug¬ 
gested with the help of wave guides by Ferraro [T2j and 
with homogeneous magnetic background fields by Dereli 
and Tucker Ca¬ 
in this paper we focus on another method for test¬ 
ing nonlinear electrodynamics and discuss its theoretical 
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foundations in detail. The basic idea is to measure the in¬ 
fluence of a (strong) background field on the propagation 
speed of light with the help of an interferometer. Such an 
experiment has been suggested and discussed already in 
five earlier papers mm- However, all of them restrict 
the theoretical discussion to the Heisenberg-Euler theory 
or, in the case of Denisov et al. [To] , to the Heisenberg- 
Euler and the Born-Infeld theory. What is still missing 
is a comprehensive derivation of the relevant equations 
that cover the whole Plebanski class. The basic idea of 
the experiment is simple. In the standard Maxwell vac¬ 
uum theory, which is linear, the superposition principle 
holds, so there is no influence of a background field on 
the propagation of light. In the nonlinear theories, how¬ 
ever, the phase velocity of light depends on the strength 
of the background field and on the propagation direction 
relative to the background field. This can be tested with 
a Michelson interferometer: If a strong background field 
is switched on and off in one interferometer arm, or if the 
whole interferometer is being rotated in a strong back¬ 
ground field, the interference pattern should change. A 
null result would place bounds on the possible deviations 
from standard Maxwell vacuum theory. 

It is the purpose of this paper to develop the theoreti¬ 
cal foundations for this experimental test for an unspec¬ 
ified nonlinear electrodynamical theory of the Plebanski 
class. We will then specify to Born, Born-Infeld and 
Heisenberg-Euler theory. 

Throughout this paper, we consider Minkowski space 
as the underlying space-time model. We work in in¬ 
ertial coordinates, so the Minkowski metric is ( rj lk ) = 
diag(l, 1,1, — 1). We use Einstein’s summation conven¬ 
tion for Latin indices taking values 1,2,3,4 and for Greek 
indices taking values 1,2,3. Indices are raised and low¬ 
ered with the Minkowski metric. We will use Gauss¬ 
ian cgs units throughout, because they are most conve¬ 
nient for our theoretical investigations. In these units, 
E, B , D and H are all measured in the same units, 
y'g/fySs). The reader can easily convert the results 
into SI units with the help of the formulas E = y/AireoEgi, 
B = v'W^Bsi, D = ^tt/VDsi and H = ^Anfj-oHsi 
m- For example, for a field X = 10 3 v /g/( v / cms) in 
Gaussian cgs units, where X = E, B, D , or U, one gets 


7 V „ 

Esi — 3 x lCfi — , B§i = 100 mT, 

m (1) 

A c A ' 

D S i = 3 x 1(T 4 H S1 = 8 x 10 4 - . 

m z m 


The paper is organized as follows. In Sec. |TT]we recall 
the basic equations for the propagation of light rays ac¬ 
cording to nonlinear electrodynamics. Then in Sec. |III 
the suggested interferometer experiment is described and 
in Sec. |IV| some particular applications are discussed. 


II. LIGHT PROPAGATION IN NONLINEAR 
ELECTRODYNAMICS 


A. The Plebanski class of nonlinear 
electrodynamical theories 

The nonlinear electrodynamical theories which are at 
the center of our examination derive from an action 

S[A m ] = ^- c j M ^ dVi . ( 2 ) 

Here j m is a given current density, A m is the electro¬ 
magnetic potential, F mn = d m A n — d n A m is the electro¬ 
magnetic field strength and C is the Larangian for the 
electromagnetic field. Then the homogeneous group of 
Maxwell’s equations is automatically satisfied, 

d [a F b c] =0. (3) 

Variation of the action with respect to the potential A m 
leads to the inhomogeneous group of Maxwell’s equa¬ 
tions, 


where 


d b H ab 



H ab 


DC 

dKb 


( 4 ) 

( 5 ) 


is the electromagnetic excitation. It is the constitutive 
law (|5j) that distinguishes different theories, while the 
Maxwell equations ([3]) and Q are always the same. 

Following Plebanski [5], we require that the electro¬ 
magnetic Lagrangian C depends on the electromagnetic 
field strength only via the Lorentz invariants 

F= l -F mn F mn and G = -±F mn F mn . (6) 


Here and in the following, the tilde denotes the Hodge 
dual, 


F ~ mn 1 ^mnab tt\ [r-r\ 

= fab- (7) 

As usual, e abcd is the totally antisymmetric Levi-Civita 
tensor with e 1234 = — 1 . Strictly speaking, only F 
is invariant under all Lorentz transformations while G 
changes sign under a parity transformation. Some au¬ 
thors restrict to Lagrangians that satisfy the equation 
C(F, G ) = C(F , —G) to assure invariance under parity 
transformations. However, for the purpose of this paper 
there is no need for this restriction. 

The Plebanski class contains, of course, the standard 
vacuum Maxwell theory which is given by the Lagrangian 

Cm{F,G) = - l -F. (8) 
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As this theory is well tested for weak fields, many authors 
restrict their work to theories where the Lagrangian satis¬ 
fies C + F/ 2 —>■ 0 for F mn —> 0. Again, for our mathemat¬ 
ical considerations there is no need for this restriction. 

For a theory of the Plebahski class the constitutive law 


([ 5 ]) can 

be written, more specifically, as 



H ab = -2 C F F ab + C G F ab 

(9) 

where 

dC dC 

(10) 


Cf= 8F ’ Ca= dG■ 


Later we will also write 


Cff 


d 2 C 
dF 2 ’ 


C-GG 


d 2 £ 
dG 2 ’ 


Cfg 


d 2 C 
dFdG ' 


( 11 ) 


Additionally wc will use a Hamiltonian formulation of 
the Plebahski electrodynamics. For the special case of 
the Born-Infeld theory, the Hamiltonian formulation can 
be found in mm Whereas the Lagrangian depends on 
the field strength, the Hamiltonian depends on the ex¬ 
citation. Quite generally, for any theory based on a La¬ 
grangian C(F mn ), the passage to the Hamiltonian for¬ 
malism can be performed whenever the constitutive law 
([ 5 ]) can be solved for F mn . The Hamiltonian is then given 
by a covariant Legendre transformation, 


U{H ab ) = ~H mn F mn - C{F ab ) (12) 


where, on the right-hand side, F ab has to be expressed 
in terms of H mn with the help of the constitutive law. 
For a theory of the Plebahski class (i.e., if the Lagrangian 
depends only on F and G) , the Hamiltonian is a function 
of the two invariants 

R=~H ab H ab and S =-^H ab H ab . (13) 


The relevant equations for the passage from the La¬ 
grangian to the Hamiltonian description are given in the 
Appendix. There we will also give a criterion that guar¬ 
antees that the constitutive law ([ 9 ]) can be solved for the 
field strength, at least locally. 


B. Three-dimensional notation of field equations 

In the following we will often use three-vector notation. 
The field strength has the three-dimensional representa¬ 
tion 

E a = F a4 , B a = ^ e a ^F M = - F“ 4 , (14) 

where e a ^ is the totally antisymmetric spatial Levi- 
Civita tensor, e 123 = 1. With E 2 = S^E^, B 2 = 
and E B = E^B^, the invariants Q read 

F=B 2 -E 2 , G = B E. (15) 


Analogously we write for the excitation 

D a = _#a4 ) H(y = 1 £aMH fh = H a4 , (16) 


which implies that the invariants (13) are given by 


R = D — H 


5 = D • H. 


(17) 


Then the constitutive law reads 

D a = -2 CpEa + CcBa , 
H a = -2£ F B a — CoE a . 


(18) 


C. Phase velocity and characteristic differential 
equation for £(F, G) theories 


The characteristic surfaces determined by a set of par¬ 
tial differential equations can be defined as the hypersur¬ 
faces along which the solutions may have discontinuities. 
As an alternative, the characteristic surfaces can also be 
defined with the help of approximate plane waves; in this 
second approach, they come about as the high-frequency 
limit of the surfaces of constant phase. In view of ap¬ 
plications to electrodynamics, the first approach is dis¬ 
cussed, e.g., in the book by Hehl and Obukhov [20]. The 
characteristic surfaces are hypersurfaces if) = constant, 
where the gradient of i/j has to satisfy, at each point of 
space-time, a fourth-order equation which is known as the 
dispersion relation or as the Fresnel equation. If viewed 
as a partial differential equation for if), this equation is 
usually called the characteristic equation or the eikonal 
equation. Using this approach, Obukhov and Rubilar pi] 
have determined the Fresnel equation (i.e., the character¬ 
istic equation) for an arbitrary C(F 7 G) theory. Earlier, 
Novello et al. [22] had found an equivalent result in a dif¬ 
ferent way. Their results show that, with the exception of 
a few special cases, theories of the Plebahski class predict 
birefringence in vacuo. For background material on bire¬ 
fringence, and bimetricity, we refer to Visser et al. [22] 
and, for the particular case of the Heisenberg-Euler the¬ 
ory, to Dittrich and Gies |2Ij and to Shore [25]. 

Here we want to briefly sketch how the Fresnel equa¬ 
tion of an arbitray C(F, G ) theory can be derived with 
the help of an approximate-plane-wave ansatz. This is 
methodically different from the work of Obukhov and 
Rubilar m and Novello et al. [22] but it leads to the 
same result. The general method goes back to Luneburg 
and is outlined, e.g., for electrodynamics in ordinary me¬ 
dia, in the book by Kline and Kay [2B]. For a discussion 
in a more general context, which includes the case to be 
considered here, we refer to Perlick [27] . 

We consider a one-parameter family of electromagnetic 
fields of the form 


F ab {x m ) = F ab (x m ) 

+ Re | {X N Fff {x m )) \ , 


N =1 


+ Re 


(19) 
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where F ab is a given background field. A is a real book¬ 
keeping parameter that is introduced in a way such that 
the high-frequency limit corresponds to A —> 0. The sum¬ 
mation sign in (191 is to be understood in the sense of 


an asymptotic series and not in the sense of a convergent 
series. While the amplitudes Fffi are in general complex, 
the eikonal function if is real. It gives the surfaces of con¬ 
stant phase, ip(x m ) = ip(x^,t) = constant. In 3-space, 
the normal to these surfaces is 


d a ip 


V(dPip) 


( 20 ) 


The phase velocity Vp can be introduced as the 3-vector 
that gives the traveling speed of such a surface in the 
direction of its normal, 


_ d a ip dip 

p (dpip)(d0ip) df [ ’ 

Feeding the ansatz into Maxwell’s equations and 
comparing equal powers of A gives a hierarchy of equa¬ 
tions. In the lowest nontrivial order, which is known as 
the geometric optics approximation, one gets a first-order 
partial differential equation for ip which is the desired 
characteristic equation. 

If this program is carried through for an £(F, G) the¬ 
ory, one finds the following result which is in agreement 
with Obukhov and Rubilar’s m- The characteristic 
equation reads 


£ F ^Mr] ij T] kl + Nrj ij F km F l m 

+ PF im F j m F kn F l n }p i p j p k p l = 0 
where pi = diip and 

M = C? F + 2 CfCfg G - - CfCgg F 
+ (F-pg ~ FffCgg ) G 2 , 

N =2£p£pF+- £ f £gg 

+ (F-fg ~ Fff£gg) F , 

P = £ff£gg — F 2 fg ■ 

If M has no zeros, ( f22j ) can be factorized as 

£ f M (afpiPj) (a^PkPe) = 0 

where 


( 22 ) 


(23) 


(24) 


of k = v ik + a A F im F k n 


for ^4 = 1,2 and 


(25) 



1 N 2 

P 

ai/2 ~ 2M ± \ 

1 AM 2 

M 


(26) 


the characteristic equation is equivalent to (which is up 
to conformal transformations in agreement with the re¬ 
sults of Novello et q/.[22]) 

a%PiPk = 0, A = 1,2 (27) 

and is sometimes called the ” light-cone condition” (com¬ 
pare for example [24jh Generalizing a standard termi¬ 
nology from electrodynamics in media, a lk and a 2 fc are 
called the optical metrics of the vacuum in the £{F , G) 
theory. If the two optical metrics do not coincide, i.e., 
if ai ^ cr 2 j there is birefringence in vacuum. If one con¬ 
siders the next order in the above-mentioned hierarchy 
of equations, one sees that the case A = 1 and the case 
A = 2 correspond to two different polarization directions. 
Note that (T\ and cr 2 are always real, because 

N 2 - AMP = (28) 

(2£p£pp — -£f£gg - PF^ +4 (^£ f £fg — PG^j 

is a sum of two squares, and that o\ and o<i depend only 
on the two field invariants F and G. In the standard 
vacuum Maxwell theory we have a± = = 0 , so these 

two functions characterize the deviation of our £(F, G) 
theory from the standard theory at the level of geometric 
optics. 

If the Lagrangian is of the special form £{F, G) = 
£(aF + f3G) with some constant factors a and /3, one has 
P = 0 and therefore ay = 0 , i.e., one polarization mode 
behaves as in the standard Maxwell vacuum theory. This 
is true, in particular, if the Lagrangian is independent of 
G. (It is also true if the Lagrangian is independent of F 
but this case was excluded by our assumption £p 0.) 
Also, it is interesting to remark that two Lagrangians £ 
and £+f3G give the same characteristic equation, i.e. the 
two cases are not distinguishable at the level of geomet¬ 
rical optics. Of course, if one restricts to parity invariant 
Lagrangians adding a term of the form j3G is forbidden. 

For some of our applications it will be desirable to write 
the optical metrics in terms of the excitation, rather than 
in terms of the field strength. It is then recommendable 
to start from a Hamiltonian formulation. It was men¬ 
tioned already at the end of Sec. |II A| tlrat the Plebanski 
class of theories can be written in terms of a Hamiltonian 
H(R,S) rather than in terms of a Lagrangian £{F,G). 
In the Appendix we derive some replacement rules of 
how the relevant Hamiltonian expressions can be found 
from the Lagrangian expressions. By applying these re¬ 
placement rules, we find that the optical metrics can be 
rewritten as 

a A = rf k + a A H im H k m (29) 

where 


In the following we restrict ourselves to Lagrangians such 
that M and £p have no zeros. This excludes some degen¬ 
erate cases which are hardly of physical interest. Then 


N 

a A = —— ± 
2 M 



P 


AM 2 M 


( 30 ) 
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with the abbreviations 


M = H ]j + 2 UrHrs S — - 'Hr'Hss R 
+ (R 2 RS ~ RrrRss) S 2 , 

N = 2 'Hr'Hrr + - Hr'Hss 

+ (R/js - H-rrH-ss ) -R, 

■P = HrrHss - Hrs ■ 


In the following it will be convenient to use the three- 
vector notation of (14) and (16). We decompose each 
of the 3-vectors E, B , D , and H into amplitude and 
direction, 


B(xt) = _B(x M , t) v(x tl , t), 
E(x >1 , t) = R(x M , t) w(x tl , t), 
H(x^,t) = H(x k ',t)r(x ,1 ,t) j 
■D(x M , t) = D(x M , t) s(x M , t), 


(32) 


where |u| = |to| = |r| = |s| = 1. The spatial and tempo¬ 
ral parts of F lrn F k rn piPki which enter into (241, can then 
be written as 


F° :m F' r n p a Pfi = B 2 p ■ p - (v ■ pf - E 2 (w ■ p) z 


F Am F^ m pp=BEp-(w x v) 

p 4 m p 4 = p 2 


(33) 


Similarly, 

H orn H 6 rrl PaPp = D 2 p p-{s pf — F[ 2 (r ■ pf 

H 4m H 0 m p 0 =DHp.(sxr) (34) 

= H 2 


Here s is a parameter along the rays which has no ob¬ 
vious physical meaning. In the following it will be con¬ 
venient to reparametrize the rays by the time coordinate 
t, cf. [23]. In order to do this, we have to assume that 
the rays of the Hamiltonian Qa are causal (i.e., timelike 
or lightlike) with respect to the Minkowski background 
metric. It was shown by Obukhov and Rubilar [2Tj that 
the optical metrics are always of Lorentzian signature, 
provided that we exclude the pathological cases where 
they degenerate. However, no convenient criterion on 
the Lagrangian C(F, G) seems to be known that guaran¬ 
tees causality of the rays with respect to the background 
metric. We will investigate this question in a separate 
paper; here we just restrict our discussion, from now on, 
to Lagrangians where the rays of the optical metrics are 
causal with respect to the background Minkowski metric. 

Then it is guaranteed that a^ 4 < 0 and we may write 
the optical metrics as 


,44 


a l APiPk = ~fr{cp4 + HJ)(cp 4: + H A ) (36) 


where 



(37) 


Equation (36) corresponds to splitting the null cone of 
the optical metric af into a future and a past cone. If 
we restrict our work here to future-oriented rays, we can 
write the characteristic equation as 


CP4 + H\ = 0 


and the canonical equations read 

dx“ d p a 

d t dp a ’ dt dx a ’ 


(38) 


(39) 


which will be used later. 


D. Ray velocity and Hamilton equations for the 
rays 

Interpreting Qa = o^APiPk as a Hamiltonian, the char¬ 
acteristic partial differential equation af dppdkip = 0 can 
be viewed as a Hamilton-Jacobi equation. The corre¬ 
sponding set of Hamilton equations, or canonical equa¬ 
tions, determines the bicharacteristic curves or rays. For 
background material on the notions of characteristics and 
bicharacteristics we refer to Courant and Hilbert [231. 

The rays are defined with respect to Q± and Q 2 sepa¬ 
rately, i.e., they depend on the polarization. The canon¬ 
ical equations read 


dx 4 

dt 


dpi 

dt 


d_H\ 
dx 4 ' 


(40) 


If a'x is known, integration of (39) gives the spatial paths 


of the rays. The first equation of (40) says that the new 


parameter t coincides with the coordinate time, while 
the second equation gives the change of the frequency of 
light. 


The ray velocity can be read from (391 


dx“ 

dt 


dH\ 

8p a 


(41) 


The phase velocity (21) can be rewritten in terms of the 
Hamiltonian as 


dx a __ OQa d p a _ 8Qa 
ds dp a ’ ds dx a 


(35) 


CP4 R 
- a P 

PaP a 


Ha 

PaP a 


P ■ 


(42) 
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Phase and ray velocity coincide if and only if 
H+ A :p0 = dK 


PaP “ " dp/3 

which is true if and only if H\ is of the form 
= f{x^,ct)y/p a p a 




(43) 


(44) 


where / is any function of the space-time coordinates. 
Equation (44) is satisfied in the usual vacuum theory of 


Maxwell but not in general in other C(F,G) theories. 
Note that (44) implies 


dx a 

d t 


d_H\ 

dp a 


f{x^,ct) p° 

\[WPf> 


(45) 


i.e., the condition Ug = Vp can hold only if dx a /dt and 
p a are parallel. 


E. Parallel electric and magnetic fields 

We consider now the special case that E and B are 
parallel, i.e., that v = w in the notation of (32). This 


case covers, of course, in particular the situation that one 
of the two field strengths, E or B , is zero. With the aid 
of the transformation (1381, described in the Appendix, 


we will then discuss, at the end of this section, the case 
that the excitations D and H are parallel. 

If we specialize (33) to the case v = w and insert the 


result into (251, the optical metrics read 


dAPiPk = -(1 ^ cr A E 2 )pl + (1 + <7 A B 2 )p 2 

- a A (B 2 + E 2 )(w ■ p) 2 . 


Hence, the Hamiltonian from (37) simplifies to 


„+ /(1 + ^aH 2 ) 

A VO- -° a e 2 ) 


[\p\ 2 - (w ■ p) 2 j + (w ■ p) 2 


and the phase velocity (42) reads 


Vp = c < 


(1 + a A B~) 
(1 ~v A E 2 ) 


1 - 


(w ■ p) 2 \ (w ■ p) 2 

~w 


I Pi 


(47) 


• (48) 


If we assume, in addition, that the unit vector of the 
background field is homogeneous, d a w^ = 0, and that 
the amplitudes of the field strengths change only in the 
direction of p, gradH oc p as well as grad A oc p , the 
canonical equations (39) reduce to 


dx“ 

dt 


H 


i f (1 + <j a B 2 ) ( ^ 

+ (w ■ p) w a } ; 


d Pa 

d t 


(49) 


OC p a ■ 


The last equation implies that the direction of p a is pre¬ 
served along the ray. 


If additionally the background fields are static, 
dE/dt = 0 and dB/dt = 0, the second equation of (40) 
reduces to 


dp4 

dt 


= 0 


(50) 


which means that, in this case, the background fields do 
not change the frequency of light. 


We are now interested in the special case that (44) 


holds which guarantees that phase velocity and ray ve¬ 
locity are equal and that dx a /dt is parallel to p a . As 
the direction of p a is preserved, the ray must then be a 
straight line. 

There are two main cases where the Hamiltonian takes 


the form of (44). First, if p 


| w, we find from (47), (48) 


and (|49|) that = c \p\, Vp = c, and dx a /dt = cp a /fpf, 


i.e., in this case the background fields have no effect. 
Second, if p ■ w = 0, one gets 


(51) 

(52) 


Vp = c 

dx a 

dt 


H y 

/ 1 + aB 2 

1 - a A E 2 

/1 + & A B 2 

1 1 - 

<7 A E 2 ’ 

/1 + o A B 2 p a 

V 1 

- a A E 2 |p| 


(53) 


This is the case which is most appropriate for the pro¬ 
posed experiment, because in this case one achieves two 
goals: the rays do not deviate from a straight line but the 
phase velocity does change in comparison to the Maxwell 
standard vacuum theory. 

Now we go over to the case that D and H are par¬ 
allel, i.e.., that r = s. With the help of (138) from the 
Appendix we find that in this case (47) and (48) have to 
be replaced with 


„ + _ /(I +d A D 2 ) 

A V {l-a A H 2 ) 


(|p| 2 - (s -p) 2 ) + (s -p) 2 (54) 


/ (1 + d A D 2 ) (, ( s-p) 2 \ ( (s p ) 2 

P V (1 - a A H 2 ) ( |p | 2 ) 


IpI 2 


|p|- 


■ (55) 


As above one gets for homogeneous and time- 
independent excitations 

dx a c f (1 + d A D 2 ) 

dt H\ { (1 - d A H 2 ) 

x [p a - (s ■ p)s a ^ + (s • p) s a | , (56) 

dp a dp4 n 

"d T*”- 

Again, the case that p is parallel to s leads to vp = c, so 
this case is of no interest for us. If, however, p ■ s = 0, 
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we get 


H+ = c\p\ 


'l + a A D 2 
1 - d A H 2 ’ 


Vp = c\ 


' 1 + a A D 2 
1 - d A H 2 


(57) 

(58) 


dx a 1 1 + d A D 2 p a 

dt c y l - <ja# 2 |p| 


(59) 


and there is no deviation of a light ray from a straight 
line. 


III. AN INTERFEROMETRIC EXPERIMENT 
FOR TESTING NONLINEAR 
ELECTRODYNAMICS 


There are two ways in which Michelson interferometry 
can be used for testing nonlinear electrodynamics. First, 
a strong background field could be applied to the light 
beam in one arm of the interferometer. One would com¬ 
pare the situation where the background field is switched 
on with the situation where it is switched off, cf. [14]. Sec¬ 
ond, one could place the whole interferometer in a strong 
background field. One would then search for changes in 
the interference pattern if the interferometer is being ro¬ 
tated. The first possibility is reasonable if one thinks of 
a large interferometer, with an arm length of several me¬ 
ters at least. The second possibility is reasonable if one 
thinks of a tabletop interferometer. As an alternative to 
using a traditional Michelson interferometer, one could 
also use a pair of optical resonators as they have been 
used for high-precision Michelson-Morley experiments in 
recent years. As these resonators have a typical size of 
only a few centimeters, one would do the experiment with 
the whole instrument placed in a background field. With 
the resonators oriented perpendicularly to each other, 
one would then compare the situation where the field 
is switched on with the situation where it is switched off, 
or one would rotate the whole instrument with keeping 
the field switched on. 

In the following we first discuss the setup of the ex¬ 
periment where a traditional Michelson interferometer is 
used and the field is placed in one arm. This is the vari¬ 
ant which brings out the basic idea of the experiment 
most clearly. Later in this section we discuss the other 
variants. 

Figure [l] shows the interferometer with the background 
field in the region denoted BF. The ray leaves the source 
S and is divided at the semipermeable mirror SPM. Af¬ 
ter reflection at the mirrors M\ and M 2 , respectively, 
both parts interfere at D. If the background field is 
switched off, both parts always travel with the stan¬ 
dard vacuum phase velocity c. If the background field 
is switched on, the part which travels along I 2 crosses 


m 2 



S Source 

D Display 

Mi Mirror 

SPM Semipermeable Mirror 
BF Background-Field 
l length of distances 

— ray 


/////////// 

D 


FIG. 1. Experimental setup 


the region BF with a different phase velocity, accord¬ 
ing to nonlinear electrodynamics. This would lead to a 
change of the interference pattern. 

We consider the background field to be static, with one 
of the four fields E, B, D , or H vanishing. Each of these 
four cases is covered by the calculations of the preceding 
section. We assume that the background field is perpen¬ 
dicular to the propagation direction of the light. We have 
seen that in this situation the ray does not deviate from 
a straight line. 


Obviously the travel times of the ray along the different 
sections are given by 

ct\ = l \, ct '2 = I 2 , ct '2 = ?2 > up tBF = Ibf ■ (60) 


Without background field the phase velocity is equal to 
c everywhere, including the region BF. The time delay 
A tj of the two arms is therefore given by 

At/ = 2(ii —t '2 t '2 — tsF) = ~(h — ^ — I'^ — Ibf) • (61) 

c 


With background field the phase velocity in the region 
BF is up which is, in general, different from c. The time 
delay At// of the two arms is therefore given by 


At// = 2(ti — t' 2 — t '2 — t B F) 

= ~(h ~ I2 ~ I2 - —Ibf) ■ 

c Up 


(62) 


The change of the interference pattern is given by the 
time difference 


At = At// — At/ 


2 Ibf 



(63) 


c 





















This leads to a line shift of 

a; At loIbf 


A = 


2 7T 


7T C 


1 - — 
Vp 


(64) 


Here w denotes the frequency of the light. Note that co 
is a constant because the background field is assumed 
static. 

We evaluate the general result for each of the four cases 
E = 0. B — (),D — 0 and H — 0. Note that in general 
E = 0 is not equivalent to D — 0 and B — 0 is not 
equivalent to H = 0. 

a) Magnetostatic field strength (E — 0): 


From (521 we find that 

v P = Csjl + a 1 / 2 B 2 

= C^1 + CTi/2 (0)^—+ 


and hence, by (641, 


A = 


uj Ibf 


1 - 


1 


7TC y yjl + <J 1/2 B 2 ) 

_ U)IbfCTi/2(0)B 2 

2 7r c + 

b) Electrostatic field strength (B — 0): 


From (521 we find that 
Vp = c 


1 


^1 - a 1/2 E 2 
= C^l+ 0-1/2 (0)^-+ 


and hence, by (641, 
w Ibf 


A = 


C7T 


(l - ^1 - a 1/2 £ 2 ) 


ujIbfcti/2^)E 2 
2 7T C 


+ ... 


c) Magnetostatic excitation (D — 0): 


From (581 we find that 


1 


Up = c 


\/l - fti/2 H 2 

= C^l + 0-1/2 (0)^-+ ... j 


and hence, by (641, 
ujIbf 


A = 


C7T 


(l-y/l-a^H 2 ) 


io Ibf&i/2(0)H 2 
2 7T C 


(65) 


( 66 ) 


(67) 


( 68 ) 


(69) 


(70) 


d) Electrostatic excitation (H — 0): 


From (58) we find that 

up = c^jl + cr 1 / 2 D 2 

= c^l + <5-1/2 (0)^-+ ...^ 


( 71 ) 


and hence, by (641, 


A _ cuZbf (i — 


^/i" 

ul BF d-i/ 2 (0)D 2 


di/2 D 2 


(72) 


2 7T C 


If one writes X for E, B, D , or H , one can combine 
all results up to first order in the form 


up =c^1 + ct 1/2 (0)^-+ ...^ 


and 


x 


A _ cuZspa- 1 / 2 (0)A 2 

2^5 h 


(73) 


(74) 


x 


Here cr 1 / 2 (0) denotes either o- 1 / 2 (0) or d- 1 / 2 (0), depend¬ 
ing on whether X is a field strength or an excitation. 
Note that 27 tc/lo is the wavelength in the case of a van¬ 
ishing background field. According to nonlinear electro¬ 
dynamics, the wavelength changes when the ray travels 
through the background field BF. This means that, if 
one substitutes the angular frequency to by the wave¬ 
length A = 2fc/uj, 


x 


A = 


Ibf& i/ 2 (0)A 2 

A 


(75) 


one has to keep in mind that A is not the wavelength of 
the light when passing through the background field but 
of the light when emitted by the source. 

The results of this section can also be applied to the 
case where the whole interferometer is inside the back¬ 
ground field. Here one does not switch on and off the 
background field but rotates the interferometer by 90° so 
that in the initial position the first arm is orthogonal to 
the field and the second arm is parallel to the field while 
in the end position it is vice versa. Then one gets instead 
of the preceding formulas the following ones: 


At = 2(/l + ?2) (1 - - 


x 


A = 


io (l i + l 2 ) a 1/2 (0)X 
2 7T C 


(76) 


(77) 
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This means that one has to replace Irf by l\ + I 2 in all 
formulas to go from the fist setup to the second one. 

As an alternative to using a traditional Michelson in¬ 
terferometer with two arms, we will now discuss a setup 
with optical resonators as it has been used frequently 
in recent years for high-precision Michelson interferome¬ 
try; see e.g. [29] and the references therein. Here one 
uses a laser which is stabilized to the eigenfrequency 
z'eigen = Nvp/ (2L) of an optical resonator, where N is 
the mode number, up is the phase velocity of light and L 
is the length of the resonator. The quality of a resonator 
is determined by its finesse F, typically F = 100 000. In 
a figurative way, a resonator may be viewed as equivalent 
to a traditional interferometer whose arm length is folded 
F times. 

If up and L undergo a change, the eigenfrequency of the 
resonator and therefore the frequency of the stabilized 
laser changes as 


8v Svp 5L 
v up L 


(78) 


If the resonator is put into a homogeneous and static E, 
B , D, or H held, with its axis perpendicular to the held, 
the phase velocity of light changes according to 


Svp 

Up 


o 1 / 2 ( 0 ) 


X 2 

~Y ' 


(79) 


if we use the approximations of (73). Asa direct measure¬ 
ment of 5v is not possible, one superimposes to the hrst 
laser a second reference laser stabilized to the eigenfre¬ 
quency t'ref of a resonator with (ideally) the same physi¬ 
cal characteristics as the hrst one. Then the difference of 
the frequencies Av := v e igen — Gel appears as the carrier 
frequency of the resulting beat. If the second resonator 
is oriented parallel to the background held and thus not 
influenced by it, this means that Av = Sv. 

For a theoretical discussion of the effect, we assume 
that L is not changed if the background held is applied. 
(Of course, for a practical realization of the experiment 
one has to take into account that the material of the 
resonator is influenced, e.g., by magnetostriction, but we 
ignore this here.) Then 


5v 

v 


(5up 

Up 


cti/ 2 (0) 


X 2 

~Y ' 


(80) 


In [29) , by averaging over many measurements it was pos¬ 
sible to determine 8v/v with an accuracy of 10 -17 . If we 
assume that the same accuracy can be reached in the 
experiment proposed here, a measurable effect requires 
that X satishes 


5vp 

Up 


£1/2(0) 


X 2 

~Y 


nr 17 . 


(81) 


In the next section we discuss the perspectives of per¬ 
forming such an experiment as a test of particular theo¬ 
ries of the Plebanski class. We compare the setup with 


the background held placed in one arm of a big Michelson 
interferometer with the setup using optical resonators. In 
the following, we refer to the hrst one as the ’’large-scale 
experiment” and to the second one as the ’’small-scale 
experiment.” 


IV. APPLICATION TO SPECIAL THEORIES 
OF THE PLEBANSKI CLASS 

A. Born-Infeld theory 

In the case of the Born-Infeld theory, Lagrangian and 
Hamiltonian are given by [2] 


/ F G 2 

C(F,G) = -b 2 0 Jl + ^--^+b 2 0 , (82) 


n(R,S) = b 2 0S Jl + ^-^-b 2 , (83) 

where bo is a new constant of Nature with the dimension 
of a held strength. The constitutive law reads 


H ab = 


jpab _ G pab 


1 


(84) 


which can be solved for the held strength, 

rmn i S 


H n 


H n 


jjmn _ 


/l + R _ si 

V b o b4 0 


(85) 


The invariants R and S are given in terms of F and G 
by 


R 

% 


bj + 4 + ^ fcg 

I 4 . £_ _ Gl 

1 ' h 2 h 4 
u 0 °0 


S = —G , 


( 86 ) 

(87) 


which implies 


1 4- b u 
1 T VI — TV 
_ 0 0 

1 + ^ 

1 + b 4 


1 4- S 

i+ bf 


I 4. R _ SI 
1 i h 2 b 4 


( 88 ) 


This leads to 


. 2 C 3 f G . 4:C 3 F G 

1~FF — - p ).— ; EfG = 


b 2 
u 0 

r 2 C F 8C 3 f G 2 
Agg = -- 


% 


b 2 o 


bl 


n, _ 2 Ur 3 S _4 Hr 3 S 

teRR 7 9 1 tlRS 


bl 


b 4 o 


R-ss = - 


2 -Hr 8R r 3 S 2 


bl 


bl 


(89) 


(90) 
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Therefore one gets for the functions <J \/2 and (Xi/ 2 , which 
give the deviation from the standard Maxwell vacuum 
theory, 


1 


CTl = a-2 = 


bl + F 


1 


CTl = (72 = — 


b 2 0 + R 


1 

+ 

°0 

1 

+ 

u 0 


(91) 


It is worth noticing that o\ = 02 holds not only in the 
Born-Infeld theory and in the standard vacuum Maxwell 
theory but also in any other theory whose Lagrangian 
differs only by a term linear in G from them. (Such 
theories, however, are often excluded because they are 
not invariant under parity transformations.) 

Additionally one can calculate the phase velocity and 
the line shift for the four static cases: 


Cases a) with E = 0 and d) with H = 0: 

Here we discuss the case of a magnetostatic field 
strength and the case of an electrostatic excitation 
together. If we use the abbreviation Y = B,D, we 
find 


Again, the limit Z —>• 0 yields us = up —> c. 

In contrast to the case above this one leads to an 
upper bound for Z. This is obvious because for 
Z —^ bo one gets us = up —> 0. In analogy to the 
other cases a background excitation slows down the 
light ray. For background fields Z > b 0 one gets an 
imaginary phase velocity, so one has to conclude 
that Z < bo . 

The line shift is given by 



(97) 


We may combine the first-order approximations of all 
preceding cases into two formulas, one for the velocity 
and one for the line shift: 

„s = „P«c(l-A) md (98) 


= ^ = (92) 

hence 


To give an example we calculate the line shift for the 
large-scale experiment for some specific values. We as¬ 
sume an accuracy of about 10 -6 line shifts. For Ibf = 
100 m and A = 1000 nm one gets: 



The limit Y —> 0 yields ug = up —> c as it has to. 

By contrast, the limit Y -» oo yields us = up —> 0, 
so one may say that the background field Y slows 
down the light ray. There is no upper bound for 
Y = B,D. 

The line shift is given by 


A = 



(94) 


Cases b) with B = 0 and c) with D = 0: 

Here we discuss the case of an electrostatic field 
strength and the case of a magnetostatic excitation 
together. If we use the abbreviation Z = E,H, we 
find 


? ' =?2= Z^’ 


(95) 


hence 


u s = Up 


Z 2 


1 - TW = C 1 - + 


Z 2 


2 bl 


(96) 


Ibf X 2 




A b 2 


-10 s 


X 2 


With these values one sees an effect if 


X > 10~ 7 b 0 


Born and Infeld conjectured that 



6 x 10 15 



(99) 


( 100 ) 


( 101 ) 


where e is the electron charge and r e is the classical elec¬ 
tron radius. Although this is only of historical interest, 
we remark that the corresponding line shift would be 


10 22 2 cms 2 

- X 2 - 

36 


( 102 ) 


If this were true we would need a field strength or an 
excitation of 


X > 6 x 10 s (103) 

Vcms 

to see an effect. For a magnetic field strength, X = B, 
this would correspond to B$i > 6 x 10 4 T in SI units, 
see ([l]). Clearly, this is not achievable in the foreseeable 
future. 

It is more interesting to see what lower bound on bo 
one could get from an experiment. Let us assume that 


X 


xlO 4 




(104) 
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which corresponds to £?si = 1T according to (JT]) . This 
is not an unrealistic value for a magnetic field to be pro¬ 
duced in a laboratory. Then a null result of ou r large- 
scale experiment would imply, according to ( 100 ), that 


bn > 1 x 10 


10 Vs 


/cms 


(105) 


For the small-scale experiment we suggest a lower field 
strength of 300 mT to prevent magnetorestriction. Then 
we find from (81) and (91) that 


b n > 7 x 10 1 


Vs 

/cms 


(106) 


which is almost 2 orders of magnitude better than the 
large-scale experiment. 


B. Born’s theory 


Born’s Langrangian |T| differs from the Born-Infeld La- 
grangian by omitting the G 2 term, 


C(F) — -blJl + ^ + b Q ■ 


(107) 


This leads to an electrodynamical theory with birefrin¬ 
gence, 


0i = 0, 


-1 


O' 2 = 


bl+F' 


(108) 


From the viewpoint of geometrical optics Born’s theory 
is a hybrid. One polarization mode behaves according to 
the standard vacuum Maxwell theory and the other one 
according to the Born-Infeld theory. This means that, if 
one filters the 02 rays out with a polarization filter, then 
one sees no difference to Maxwell, and if one filters the 
u 1 rays out, then one sees no difference to Born-Infeld. 
As a consequence, the results of Sec. IV A are also valid 
for the cr 2 rays in Born’s theory. 


C. Series expansions for electrodynamics with 
arbitrary Lagrangian 


If we are interested only in first-order deviations from 
Maxwell’s theory, we may express the Langrangian in 
terms of a series expansion with respect to F/A and G/A 
up to second order, where A is a constant with the dimen¬ 
sion of a held strength squared. Introducing A is neces¬ 
sary because only for dimensionless terms is it meaningful 
to say that they are small without referring to a partic¬ 
ular system of units. In the Born-Infeld theory, e.g., we 
choose A = 6q. 

The series expansion of the Lagrangian reads 


C = a + 


F 

A 


+ P 

2 


G 

A 


+ 


+ 7 
1 



FG 


■7 

3 



(109) 


Note that /? and 7 are zero if the theory is invariant under 
2 2 

parity transformations. One can assume the validity of 
the following bookkeeping system for the smallness of 
terms, where ~ means that terms are of the same order. 


• a ” 

G. 


7 • • • as well as F mn ~ F rnn , p 


F/A and G/A are dimensionless with F/A ~ G/A. 

From the first ord er of R = R(F,G ) and S = 
S(F,G) [cf. (139) to (145)] one gets fiF/A ~ 

AR/p. 


With the help of (139) to (145)) one can now calculate 
R and S as series in F and G. Additionally one can then 
calculate the inverted series, i.e. F and G as a series 
in R and S. This step allows us then to calculate the 
Hamiltonian as a function of R and S. The result of this 
calculation is 


U = ~a + B I l+B ! i 
1 A 2 A 

' td\ 2 

■<7| 5 


A 


+ C™ 

2 ,42 


C 

3 


T +■ 


with the following coefficients: 

A 


A” 1 := 


/r + 4/5 


B~-P\ B:=i 3 ; 
1 1 2 2 


( 110 ) 


C := —p 7 + 2/3/r 7 - 4/V/3 7 + 8/? 2 /3 2 7 - 


2 1 


12 2 


12 3 


12 1 


C:= 

2 


-8/? 3 /?7-16/3S) / (V+4/T 

12 2 1 1 / / \2 1 


—/5 4 7 + 4/5/5 3 7 — 16/?/5 3 7 + 24/5 2 /9 2 7— 


12 3 


(111) 


-16/3 3 /37 + 64/3 3 /37- 16/3 4 7 ) / (V+4/3' 

123 121 1 2/' \2 1 

C := ( —/3 4 7 - 8/5/5 3 7 + 8/3 2 /3 2 7 - 64/3 2 /3 2 7 + 


2 3 


122 


12 3 


12 1 


+32/? 3 /3 7 - 16/5 4 7 ) / (/3 2 + 4/5' 

122 1 3/ / \2 1 


This leads to some additions to the bookkeeping system: 

• From pF/A ~ AR/P one gets F/A ~ R/A as well 

i i 

as G/A ~ S/A. 

• It is easy to see that 

i i i i 

holds. 


As a result of this bookkeeping system one sees that, if 
one can neglect in the Lagrangian terms of a certain order 
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in F/A and G/A, then one can neglect in the Hamiltonian 
terms of the same order in R/A and S/A. 

For the case of a parity-invariant Lagrangian, 
C(F, G) = C{F, — G), the coefficients of the Hamiltonian 
become very simple: 


B = ~P\ 

1 1 

? = -7 ! 


B = 
2 

c = 

2 


~P = 0; 

2 

-7 = 0 ; 

2 


c 

3 


"7, 

3 


( 112 ) 


A = 


4/3 2 

i 

X 


If, in addition, the first-order approximation of the La¬ 
grangian coincides with the standard Maxwell one (| 8 j) — 
which is true if /3 = —A /2 — one gets 


Magnetostatic field strength case (E = 0): 


Vp — c + c 


W7 7 


2 l A/3 2 A/3 

\ l i 




\ 


I67 2 + 47 2 — 877 + 7 2 

1 _2_ 13 3 

X/ 3 2 

1 


1 R 2 / 2 7 7 

a= 1bfB x+ — 

A l A/3 2 A/3 

\ 1 1 


±- 


! \ 


I67 2 + 47 2 — 877 + 7 2 

1 _2_ 13 3 

A 2 /3 2 

1 


(116) 


+ •• 


(117) 


A = A. 


(113) Electrostatic field strength case (B = 0): 


This is, in particular, the case for the theories of Born, 
Born-Infeld, and Heisenberg-Euler. 

From the Lagrangian (or the Hamiltonian, respec¬ 
tively) one gets the “deviation coefficients” 07/ 2 in the 
zeroth order of approximation with respect to F and G 
(or R and S , respectively): 


up = c + c 


£ 2 /27 7 






+ 


2 \ A/3 2A(3 

\ 1 1 


\ 


I67 2 + 47 2 — 877 + 7 2 

1 _2_ 13 3 

A 2 P 2 
1 


(118) 


+ •• 


2 7 


7 

3 


ai/2{0) -Jp + 2A/3 
1 1 


1 

*5 


\ 


I67 2 + 47 2 — 877 + 7 2 

1 2 13 3 


A 2 (3 2 
1 


2 C C 

" ,/2<0) = Is + 23b 
1 1 


1 

zb — 
2 


\ 


16C 2 + AC 2 - 8 CC + C 2 

1 _2_ 13 3 

A 2 B 2 


(114) 


■ (H5) 


Obviously the zeroth order approximation of 07/2 gives 
the first-order approximation of the optical metric for the 
deviation from Maxwell’s theory and therefore also of the 
phase velocity. 

The Born-Infeld theory, e.g., yields ay(0) = < 72 ( 0 ) = 
cti(0) = <j 2 (0) = —-1 /fog > so one recovers the values cal¬ 
culated above. Additionally one sees that the approxi¬ 
mation procedure does not destroy the absence of bire- 
frigence in the given order of approximation. 

O ne g ets the res ults f or th e fou r cases described in 
Sec. [HI if one feeds (|115|) and (|114|) into ([73|) and (J74[) : 


A = 


Ibf E 


2 / 2 7 7 


A \AP 2Ap 
1 1 


±1 


A 


I 67 2 + 47 2 — 877 + 7 2 

1 _2_ 13 3 

A 2 / 3 2 
1 


+ 


Magnetostatic excitation case (D = 0): 


Up = c + c 


±- 


tt2 / 2C C 

n \ 1,3 


2 \ AB 2AB 

1 1 


A 


16 C 1 + 4 C 2 - 8 CC + C 2 

1 2 13 3 


A 2 B 2 


A = 


, TT2 (2C C 

Iff n / i 


A \ AB 2AB 

\ 1 1 


zb- 


A 


16 C 2 + 4 C 2 - 8 CC + C 2 

1 _2_ 13 3 

A 2 b 2 


(119) 


( 120 ) 


( 121 ) 
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Electrostatic excitation case (H = 0): 


Vp =c + c 


D 2 


± : 


: \ 


( 2C 

(4 


+ 


C 

_3_ 

2 AB 


16G 2 + 4C 2 - 8CC + G 2 

1 _2_ 13 3 

A 2 b 2 

1 


( 122 ) 


A = 



A 


16G 2 + 4G 2 - 8 CC + C 2 

1 _2_ 13 3 

A 2 b 2 


(123) 


In principle it is easy to obtain further orders of ap¬ 
proximation, but the resulting terms are expected to be 
very small and will not be worked out here. 


For the four possibilities for the background held de¬ 
scribed in Sec. |III| we get, using again the abbreviation 


X = E,D,B,H, 


Vp(<7i) = C |^1 

vp(0- 2 ) = c ( 1 


7A X 2 \ 

“eT ) 


4A X 2 \ 

) 


A(cri) 


A(<7 2 ) 


UAIbfX 2 

~Tb 2 ~ + 


8 A IbfX 2 

A £ 0 2 + 


(130) 


Using the same setup as before for the large-scale ex¬ 
periment, with Ibf = 100 m and A = 1000 nm one gets 


A(cti) 

A(ct 2 ) 


-2 x 10~ 23 X 2 


cms 2 

g 


-1 x 10~ 23 X 2 


cms 2 

g 


(131) 


D. The Heisenberg-Euler theory 


Here we give an example for the procedure described 
in the preceding section. For small values of the held 
strength the Heisenberg-Euler theory can be described 
by the following Lagrangian Eiim which results from a 
series expansion with respect to F and G: 


where 


1 F 


C = EU---^+ Al—r+7—r 


2 E 2 


A = 


F 2 


G 2 


F 4 


E 4 


he 


907re 2 


= 0.7363 


En = 


2 4 

me 


= 6.048 x 10 15 


Vs 

/cms 


(124) 


(125) 

(126) 


Here e is the electron charge, m is the electron mass, c is 
the speed of light and h is Planck’s constant. 

So the coefficients in (109) are 


a = 0, P = -B = - E ^-, P = 0, 

1 1 Z 2 

1 = = A E%, 7 = -G = 0, (127) 

1 1 2 2 

7 =-G = 7AE 2 . 

3 3 

and 

A = A = E 2 . (128) 


Hence 


o-r(0) = di(0) = --^2“. 

o a 

o - 2 ( 0 ) = ct 2 ( 0 ) = -^2 • 

En 


(129) 


Therefore one needs a held strength or an excitation of 

X > 3 x 10 8 A (132) 

Veins 

to see any effect. This is clearly not achievable with 
present or near-future instruments. 

A similar calculation shows that for the small-scale 
experiment a held about 2 orders of magnitude smaller 
would be sufficient. However, even in this case one would 
need a held of more than 10 6 =100 T to see an effect. 


V. CONCLUSIONS AND DISCUSSION 

Since Born and Infeld created their “new held theory” 
of electromagnetism [2], different nonlinear modifications 
of vacuum electrodynamics on the basis of a Lagrangian 
C(F, G ) have been discussed, where usually one considers 
only those theories that reproduce the standard vacuum 
Maxwell theory in sufficiently weak helds. All these new 
electrodynamical theories have in common that they pre¬ 
dict that light travels along the null cones of two optical 
metrics, one for each polarization state, where at least 
one of them differs from the vacuum Maxwell light-cone. 
At the same time they introduce at least one new dimen- 
sionfull constant of Nature. 

While in the standard vacuum Maxwell theory the su¬ 
perposition principle holds, this is no longer true in other 
C(F, G ) theories. As a consequence, an electromagnetic 
background field would have an effect on the propagation 
of electromagnetic waves and thus, in particular, on the 
phase velocity of light. This is reflected by the fact that 
the optical metrics depend on the background field. The 
best technique for measuring small changes in the phase 
velocity of light with high accuracy is interferometry. In 
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this paper we worked out the mathematical details for 
using interferometry as a test of C(F, G ) theories. 

In cases where the constants of Nature that enter into 
the theory are known, as e.g. in the Heisenberg-Euler 
theory, an interferometric experiment could be used for 
confirming the theory by verifying the prediction. If in¬ 
stead the constants of Nature that enter into the the¬ 
ory are not known, as e.g. in the Born-Infeld theory, a 
null result of the experiment would give bounds on these 
constants. Our estimates demonstrate that, with realis¬ 
tic (magnetic) fields, an interferometric experiment could 
place significant bounds on the Born-Infeld constant bo. 

Unfortunately, in the case of the Heisenberg-Euler 
theory our estimates seem to indicate that a confirma¬ 
tion of the theory is not realizable with electromag¬ 
netic fields that can be achieved in present-day exper¬ 
iments. However, it might be possible to considerably 
enhance the sensitivity by using time-dependent back¬ 
ground fields, rather than the static fields we have con¬ 
sidered for our numerical estimates. For the case of test¬ 
ing the Heisenberg-Euler theory with an interferometer 
of the size of a gravitational wave detector, this possibil¬ 
ity was discussed in detail recently by Grote [15] . The 
idea is to change the background field periodically with 
a frequency ui, e.g. by rotating a permanent magnet. 
As long as u> is small in comparison to the frequency of 
the laser light used in the interferometer, our equations 
could still be used for this situation in the sense of an 
adiabatic approximation. If the laser light is polarized, 
rotating the background field would lead to a periodi¬ 
cally varying signal according to any theory that predicts 
birefringence in vacuo. (Unfortunately, this excludes the 
Born-Infeld theory.) By choosing long integration times 
- Grote suggests to run the experiment for a year - 
one could improve the statistics in such a way that it 
might be possible to reach the sensitivity for testing the 
Heisenberg-Euler theory. A similar analysis has not been 
carried through for the small-scale experiment so far. We 
will leave this for other authors, as it goes beyond the 
scope of the present paper which was to lay the theoret¬ 
ical foundations of the experiment in the context of an 
arbitrary £(F, G ) theory. 

Finally, we add a remark on pulsed background fields. 
Pulsed magnetic fields and also laser pulses (pulsed null 
fields) can be produced with considerably higher field 
strengths than static or slowly varying fields. For exam¬ 
ple, pulsed magnetic fields of ~ 100 T have already been 
produced in the laboratory. However, these fields per¬ 
sist only for short times, so the adiabatic approximation 
would not be valid which makes the theory considerably 
more difficult. Moreover, there are several technical ob¬ 
stacles. For example, we see major experimental difficul¬ 
ties towards a realization of the small-scale experiment 
with (pulsed) magnetic fields of ~ 100 T because of mag¬ 
netostriction. Also, for the experiment with a pulsed null 
field as a background one would wish to have the pulse 
traveling in the same direction as the laser beam in the 
interferometer, to make sure that the latter does not de¬ 


viate from a straight line. This cannot be done without 
changing the geometry of the interferometer, neither for 
the small-scale nor for the large-scale experiment. For 
these reasons, we have restricted our specific calculations 
to time-independent background fields (which includes 
the case of slowly varying fields in the sense of an adia¬ 
batic approximation). 
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APPENDIX: HAMILTONIAN FORMALISM IN 
TERMS OF THE EXCITATION 


First we give a necessary and sufficient condition for 
the constitutive law § to be locally solvable for F ab . By 
the implicit function theorem, this is true if the Jacobian 
of the map from the field strength 6-vector to the exci¬ 
tation 6-vector is nonzero After dividing by the factor 
(4 C 2 f + £q) 2 , which is nonzero unless the Lagrangian 
is constant and thus trivial, we find that this condition 
reads 


4£| + _ 4 (£ 2 + 4G 2 ) c ff C gg + 4 ( F 2 + 4 G 2 ) C? FG 

+8F£ f £ff + 16 G£f£-fg + 4 FCgg^g ~ “^FCpCcG 
— 8G£ff£g + 2G£g£gg 7^ 0 • (133) 


It is easy to see that this condition is satisfied, for all 
field configurations, in the Born theory and also in the 
Born-Infeld theory. For the Heisenberg-Euler Lagrangian 
(1241 it is true as well, where we have to observe that this 
second-order theory is valid only as long as the magnitude 
of the field strength is small in comparison to Eq. 

Whenever the constitutive law ([5]) can be solved for 
F rnn . we can pass to a Hamiltonian description by a Leg¬ 
endre transformation (12). In this appendix we derive 
some relevant equations of the Hamiltonian formalism 
that will be used in the body of the paper, based on an 
analogue formalism that was developed already by Born 
and Infeld [2] for their special theory. 

From © and ([5]) we find 


&H 
d Hg 


— —F- 

— * ij 


( 134 ) 
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which is the Hamiltonian version of the constitutive law. 
In the case of vanishing sources, j m = 0, the Maxwell 
equations read 

d n H mn = 0 and d [a F bc] = 0. (135) 

These two equations can be equivalently rewritten as 

d [a H bc] = 0 and d n F mn = 0. (136) 


Comparison of ([5]) and (135) on one side and (134) and 
(1361 on the other side demonstrates that the source-free 


theory is invariant under a duality rotation 

pmn ^ pmn £^H. 


(137) 


In 3-vector notation, F mn means E a c —>• H a and 

B a <->■ -D a . Clearly, F mn -A H mn implies 


pmn ^ _pn 


F^R, S. 


(138) 


If we start from the Lagrangian £(F mn ) and work out 
all relevant equations of the theory in terms of the field 
strength, we get the relevant equations in terms of the 
excitation simply by applying the replacements (137) and 
( 138| ). Note that this method works only in the case 
of vanishing sources, j m = 0, but for any Lagrangian 
£{F mn ) for which the constitutive law ([5]) can be solved 
for F rnn . 

We now specify to a Lagrangian of the Plebanski class. 
We recall that in this case the constitutive law reads 


H ab = -2 C F F ab 


Similarly, (134) specifies to 

F a b = 2 Hr H ab - Hs H ab . 


Inserting (139) into (12) yields 


H(R, S ) = 2 £ f F + 2 C g G - £(F, G) 


while inserting (140) into (12) yields 


H{R, S ) = 2 HrR + 2 H s S - £(F, G ). 
From these two equations we read that 

C f F + C g G = HrR + H s S . 
Also, from ( 139| ) we find immediately that 

R = (—4 C? F + C 1 2 3 4 5 6 7 8 9 10 11 g ) 2 F - 8 C f C g G , 
S = (—A£ 2 f + C G ) 2 G + 2 C f C g F . 

Similarly, from (|140[) we find that 


F = (-4 H 2 r + H 2 s ) 2 R- 8HrH s S , 
G = (-4 H\ + H 2 s ) 2 S + 2 HrH s R . 


(140) 


(141) 


(142) 


(143) 


(144) 


(145) 


CgF 


ab 


(139) 


In Sec. IV C| the equations ( |139[ ) to (145) are used for 
calculating series expansions of the Lagrangian and the 
Hamiltonian theory up to second order in F and G. This 
enables one to calculate the first post-Maxwellian results 
of the discussed experiment for an arbitrary Lagrangian 
of the Plebanski class. 
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